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Abstract 

Evidence has recently emerged for a hidden symmetry of scattering amplitudes in 
M = 4 super Yang-Mills theory called dual conformal symmetry. At weak coupling the 
presence of this symmetry has been observed through five loops, while at strong coupling 
the symmetry has been shown to have a natural interpretation in terms of a T-dualized 
AdS*,. In this paper we study dual conformally invariant off-shell four-point Feynman 
diagrams. We classify all such diagrams through four loops and evaluate 10 new off-shell 
integrals in terms of Mellin-Barnes representations, also finding explicit expressions for 
their infrared singularities. 
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1. Introduction 

Recent work on M = 4 super Yang-Mills theory has unlocked rich hidden structure 
in planar scattering amplitudes which indicates the exciting possibility of obtaining exact 
formulas for certain amplitudes. At weak coupling it has been observed at two and three 
loops fli~4] , |r7|| that the planar four-particle amplitude satisfies certain iterative relations 



which, if true to all loops, suggest that the full planar amplitude A sums to the simple 
form 

\og(A/A t ree) = (IR divergent terms) + log 2 (t/s) + c(A) H , (1.1) 



where /(A) is the cusp anomalous dimension, s and t are the usual Mandelstam invariants, 
and the dots indicate terms which vanish as the infrared regulator is removed. Evidence 
for similar structure in the five-particle amplitude has been found at two loops [pl| , |22|| . At 



strong coupling, Alday and Maldacena |S7| have recently given a prescription for calculating 
gluon scattering amplitudes via AdS/CFT and demonstrated that the structure ( |TTT| ) holds 
for large A as well. 

An important role in this story is evidently played by a somewhat mysterious symme- 



try of M = 4 Yang-Mills theory which has been called 'dual conformal' symmetry in [29]. 
This symmetry, which is apparently unrelated to the conventional conformal symmetry 
of M = 4 Yang-Mills, acts as conformal transformations on the variables Xi = fcj — fcj+i, 
where fcj are the cyclically ordered momenta of the particles participating in a scattering 
process. It is important to emphasize that dual conformal invariance is a property of planar 
amplitudes only. Although somewhat mysterious at weak coupling, the Alday-Maldacena 
prescription provides a geometrical interpretation which makes dual conformal symmetry 
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manifest at strong coupling. One of the steps in their construction involves T-dualizing 
along the four directions of AdS$ parallel to the boundary, and dual conformal symmetry 
is the just isometry of this T-dualized AdS$. 

The generalized unitarity methods P,fF|,|9l ]l6| , [18| which are used to construct the 
dimensionally-regulated multiloop four-particle amplitude in M = 4 Yang-Mills theory 
express the ratio A/ At re e as a sum of certain scalar Feynman integrals — the same kinds of 
integrals that would appear in <p n theory, but with additional scalar factors in the numera- 
tor. However dimensional regularization explicitly breaks dual conformal symmetry, so the 



authors of |?9| used an alternative regularization which consists of letting the four external 
momenta in these scalar integrals satisfy kf — fi 2 instead of zero. Then they observed that 
the particular Feynman integrals which contribute to the dimensionally-regulated ampli- 
tude are precisely those integrals which, if taken off-shell, are finite and dual conformally 
invariant in four dimensions (at least through five loops, which is as far as the contributing 
integrals are currently known [|T|,|8|, p!7| , p4| , |26[| ) . 

Off-shell dual conformally invariant integrals have a number of properties which make 
them vastly simpler to study than their on-shell cousins. First of all they are finite in four 
dimensions, whereas an L-loop on-shell dimensionally regulated integral has a complicated 
set of infrared poles starting at order e~ 2L . Moreover in our experience it has always proven 
possible to obtain a one-term Mellin-Barnes representation for any off-shell integral, several 
examples of which are shown explicitly in section 4. In contrast, on-shell integrals typically 
can only be written as a sum of many (at four loops, thousands or even tens of thousands 
of) separate terms near e = 0. It would be inconceivable to include a full expression for 
any such integral in a paper. 

Secondly, the relative simplicity of off-shell integrals is such that a simple analytic 
expression for the off-shell L-loop ladder diagram was already obtained several years 
ago [§,§J^] (and generalized to arbitrary dimensions in |13|]). For the on-shell ladder dia- 
gram an analytic expression is well-known at one-loop, but it is again difficult to imagine 
that a simple analytic formula for any L might even be possible. 

Finally, various off-shell integrals have been observed to satisfy apparently highly 
nontrivial relations called 'magic identities' in (23|| . There it was proven that the three- 
loop ladder diagram and the three-loop tennis court diagram are precisely equal to each 
other in four dimensions when taken off-shell. Moreover a simple diagrammatic argument 
was given which allows one to relate various classes of integrals to each other at any number 
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of loops. No trace of this structure is evident when the same integrals are taken on-shell 
in 4 — 2e dimensions. 

Hopefully these last few paragraphs serve to explain our enthusiasm for off-shell in- 



tegrals. Compared to our recent experience P5|j3lH with on-shell integrals, which required 
significant supercomputer time to evaluate, we find that the off-shell integrals we study 
here are essentially trivial to evaluate. 

Unfortunately however there is a very significant drawback to working off-shell, which 
is that although we know (through five loops) which scalar Feynman integrals contribute to 
the dimensionally-regulated on-shell amplitude, we do not know which integrals contribute 
to the off-shell amplitude. In fact it is not even clear that one can in general provide a 
meaningful definition of the 'off-shell amplitude.' Taking kf = fx 2 in a scalar integral 
seems to be a relatively innocuous step but we must remember that although they are 
expressed in terms of scalar integrals, the amplitudes we are interested in are really those 
of non-abelian gauge bosons. In this light relaxing the on-shell condition kf = does not 
seem so innocent. If it is possible to consistently define a general off-shell amplitude in 
M = 4 Yang-Mills then we would expect to see as fi 2 — > the universal leading infrared 
singularity 

log (A/Atree) = - — log (// / st) + less singular terms, (1-2) 



where /(A) is the same cusp anomalous dimension appearing in If an equation of 

the form (|1 . 2| ) could be made to work with off-shell integrals, it would provide a method 
for computing the cusp anomalous dimension that is vastly simpler than that of reading it 



off from on-shell integrals as in [24,25,31 



The outline of this paper is as follows. In section 2 we review the definition and dia- 
grammatic properties of dual conformal integrals. In section 3 we present the classification 
of off-shell dual conformal diagrams through four loops. In section 4 we evaluate 10 new 
dual conformal integrals in terms of Mellin-Barnes representations (finding two new 'magic 
identities') and present explicit formulas for their behaviour in the infrared limit [x 2 — > 0. 



2. Properties of Dual Conformal Integrals 

We begin by reviewing the definition and diagrammatic properties of dual conformal 



integrals following |23|,£9|. By way of illustration we consider first the one-loop diagram 
shown in Fig. 1. Black lines depict the underlying scalar Feynman diagram, with each 
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internal line associated to a 1/p 2 propagator as usual. Each dotted red line indicates a 
numerator factor (pi + p 2 + ■ ■ ■ + p n ) 2 where the pi are the momenta flowing through the 
black lines that it crosses. Of course momentum conservation at each vertex guarantees 
that a numerator factor only depends on where the dotted red line begins and ends, not 
on the particular path that it traverses through the diagram. 



Fig. 1: The one- loop scalar box diagram with conformal numerator factors 
indicated by the dotted red lines. 



We adopt a standard convention (see for example ||14|| ) that each four-dimensional 
loop momentum integral comes with a normalization factor of 1/iir 2 . Thus the diagram 
shown in Fig. 1 corresponds to the integral 



l {1) (h,k 2 ,k 3 ,k 4 ) = J 



d 4 pi 



(h + k 2 ) 2 (k 2 + k 3 



2 Pi(pi ~ h) 2 (pi - h - k 2 ) 2 (pi + k 4 ) 2 ' 



VK 



(2.1) 



We regulate this infrared divergent integral by taking the external legs off-shell, choosing 
for simplicity all of the 'masses' kf = [i 2 to be the same. A different possible infrared 
regulator that one might consider would be to replace the 1/p 2 propagators by massive 



propagators 1 / (p z 



m 



but we keep all internal lines strictly massless. 



Following p3| , |29[1 we then pass to dual coordinates Xi by taking 



k\= x 



12; 



k 2 — x 2 3, k% — £34, k 4 — X41, pi — £15, 



(2.2) 



where Xij = Xi — Xj, so that Q2.1|) becomes 

l ( - 1) (x ll x 2l x 3l x 4 ) = J 



24 



• 9 2 2 2 2' 

tn x 15 x 2 5X3 5 x 4 5 



(2.3) 



This expression is now easily seen to be invariant under arbitrary conformal transforma- 



tions on the Xi. This invariance is referred to as 'dual conformal' symmetry in |23[ because 
it should not be confused with the familiar conformal symmetry of M = 4 Yang-Mills. The 
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Fig. 2: The one- loop scalar box with dotted red lines indicating numerator 
factors and thick blue lines showing the dual diagram. 

coordinates x% here are momentum variables and are not simply related to the position 
space variables on which the usual conformal symmetry acts. 

To analyze the dual conformal invariance of a diagram it is convenient to consider its 
dual diagramR In Fig. 2 we have labelled the vertices of the dual diagram to Fig. 1 in 



accord with the expression ( |2.3| ). The numerator factors correspond to dotted red lines as 
before, while denominator factors in the dual expression ( |2.3| ) correspond to the thick blue 
lines in the dual diagram. Neighboring external faces are not connected by thick blue lines 
because the associated propagator is absent. 

With this notation set up it is easy to formulate a diagrammatic rule for determining 



whether a diagram can correspond to a dual conformal integral |2J]. We associate to each 
face in the diagram (i.e., each vertex in the dual diagram) a weight which is equal to the 
number of thick blue lines attached to that face minus the number of dotted red lines. 
Then a diagram is dual conformal if the weight of each of the four external faces is zero 
and the weight of each internal face is four (to cancel the weight of the corresponding 
loop integration f d 4 x). Consequently, no tadpoles, bubbles or triangles are allowed in 
the Feynman diagram, each square must be associated with no numerator factors, each 
pentagon must be associated with one numerator factor, etc. 

We distinguish slightly between dual conformal diagrams and dual conformal integrals. 
The latter are all diagrams satisfying the diagrammatic rule given above. However in |29| 
it was pointed out that not all dual conformal diagrams give rise to integrals that are finite 
off-shell in four dimensions. Those that are not finite can only be defined with a regulator 



1 The fact that nonplanar graphs do not have duals is consistent with the observation that 
dual conformal symmetry is apparently a property only of the planar limit. 



(3)b (3)c 

Fig. 3: Two examples of three- loop dual conformal diagrams. 

(such as dimensional regularization) that breaks the dual conformal symmetry and hence 
cannot be considered dual conformal integrals. 

In Fig. 3 we show two diagrams that are easily seen to be dual conformal according 
to the above rules. The diagram on the left is the well-known three- loop tennis court 
j(3)6 Thedia gram on the right demonstrates a new feature that is possible only when 



the external lines are taken off-shell. The dotted red line connecting the top external face 
to the external face on the right crosses only one external line and is therefore associated 
with a numerator factor of kf = fx 2 . Such an integral is absent when we work on-shell. Of 
course absent does not necessarily mean that an integral vanishes if we first calculate it 
for finite /i 2 and then take y? — > 0. Indeed we will see below that ~ ln 3 ^ 2 ) in the 
infrared limit. 

An important and well-known feature of four-point dual conformal integrals is that 
they are constrained by the symmetry to be a function only of the conformally invariant 
cross-ratios 

2 2 2 2 

M= ^§i, ^ = ^23. (2.4) 

x 13 x 24 x 13 x 24 

Since we have chosen to take all external momenta to have the same value of kf = fi 2 , we 
see that u and v are are actually both equal to 

st 

where s = (pi + P2) 2 and t = {j>2 + P3) 2 are the usual Mandelstam invariants. Therefore 
we can express any dual conformal integral as a function of the single variable x. One 
immediate consequence of this observation is that any dual conformal integral is invariant 
under rotations and reflections of the corresponding diagram, since x itself is invariant 
under such permutations. A second consequence is that any degenerate integral (by which 
we mean one where two or more of the external momenta enter the diagram at the same 
vertex) must evaluate to a constant, independent of x. 
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3. Classification of Dual Conformal Diagrams 

3.1. Algorithm 

Let us now explain a systematic algorithm to enumerate all possible off-shell dual 
conformal diagrams. We use the graph generating program qgraf to generate all scalar 
ipi four-point topologies with no tadpoles or bubbles, and throw away any that are non- 
planar or contain triangles since these cannot be made dual conformal. After these cuts 
there remain (1, 1, 4, 25) distinct topologies at (1, 2, 3, 4) loops respectively respectively. We 
adopt the naming conventions from reference [J24] which displays 24 out of the 25 four-loop 
topologies, omitting the one we call h in Fig. 9 since it vanishes on-shell in dimensional 
regularization. Note also that the topologies es and c\ shown there are actually the same. 

The next step is to try adding numerator factors to render each diagram dual confor- 
mal. Through three loops this is possible in a unique way for each topology, but at four 
loops there are three topologies (shown in Fig. 4) that cannot be made dual conformal at 
all while six topologies (&i, c, <i, e, e2 and /, shown below) each admit two distinct choices 
of numerator factors making the diagram dual conformal. 




Fig. 4: These are the three planar four-point IPI four-loop tadpole-, bubble- 
and triangle-free topologies that cannot be made into dual conformal diagrams 
by the addition of any numerator factors. In each case the obstruction is that 
there is a single pentagon whose excess weight cannot be cancelled by any 
numerator factor because the pentagon borders on all of the external faces. 
(There are no examples of this below four loops.) 

We remark that we exclude from our classification certain 'trivial' diagrams that can 
be related to others by rearranging numerator factors connected only to external faces. For 
example, consider the two diagrams shown in Fig. 5. Clearly both are dual conformal, but 
they differ from each other only by an overall factor of x = fi^/st. In cases such as this we 
include in our classification only the diagram with the fewest number of \j? powers in the 

2 We do not know of a general proof that no one-particle reducible diagram can be dual 
conformal but it is easy to check through four loops that there are no such examples. 
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Fig. 5: Two dual conformal diagrams that differ only by an overall factor. 
As explained in the text we resolve such ambiguities by choosing the integral 
that is most singular in the fi 2 — > limit, in this example eliminating the 
diagram on the right. 

numerator, thereby choosing the integral that is most singular in the \j? —>■ limit. In the 
example of Fig. 5 we therefore exclude the diagram on the right, which actually vanishes 
in the fi 2 — > limit, in favor of the diagram on the left, which behaves like ln 2 (/x 2 ). 

Another possible algorithm, which we have used to double check our results, is to first 
use the results of ||12|| to generate all planar 1PI vacuum graphs and then enumerate all 
possible ways of attaching four external legs so that no triangles or bubbles remain. 



3.2. Results 

Applying the algorithm just described, we find a total of (1,1,4,28) distinct dual 
conformal diagrams respectively at (1,2,3,4) loops. While (1,1,2,10) of these diagrams 
have appeared previously in the literature on dual conformal integrals [|23|,[24],^B| , the re- 
maining (0, 0, 2, 18) that only exist off-shell are new to this paper. We classify all of these 
diagrams into four groups, according to whether or not they are finite in four dimensions, 
and according to whether or not the numerator contains any explicit factors of /i 2 . Hence 
we define: 

Type I diagrams are finite in four dimensions and have no fi 2 factors. 

Type II diagrams are divergent in four dimensions and have no [i 2 factors. 

Type III diagrams are finite in four dimensions and have [i 2 factors. 

Type IV diagrams are divergent in four dimensions and have // 2 factors. 

Type I and II diagrams through five loops have been classified, and some of their 
properties studied, in [23,24,26 ,^9|. In particular, it has been observed in these references 
that it is precisely the type I integrals that contribute to the dimensionally regulated 
on-shell four-particle amplitude (at least through five loops). We display these diagrams 
through four loops in Figs. 6 and 7. The new type III and IV diagrams that only exist 
off-shell are shown respectively in Figs. 8 and 9. Each diagram is given a name of the 
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form X^ 1 where L denotes the number of loops and i is a label. The one- and two-loop 
diagrams are unique and do not require a label. Below we will also use Z( L ) to refer to the 
L-loop ladder diagram (specifically, ZW, 1^ 2 \ and for L = 1,2, 3, 4). 

We summarize the results of our classification in the following table showing the 
number of dual conformal diagrams of each type at each loop order: 



L 


I 


II 


III 


IV 


1 


1 











2 


1 











3 


2 





2 





4 


8 


2 


9 


9 



(3.1) 
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Fig. 6: Type I: Here we show all dual conformal diagrams through four 
loops that are finite off-shell in four dimensions and have no explicit nu- 
merator factors of fi 2 . These are precisely the integrals which contribute to 
the dimensionally-regulated on-shell four-particle amplitude p||^, p~7| , |2"4}j . For 
clarity we suppress an overall factor of st in each diagram. 
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"V 



d5 



Fig. 7: Type II: These two diagrams have no explicit factors of \j? in the 
numerator and satisfy the diagrammatic criteria of dual conformality, but the 
corresponding off-shell integrals diverge in four dimensions ||29|1 . (There are 
no examples of this type below four loops). 







\ e2 









V d 


b2 V el 



Fig. 8: Type III: Here we show all diagrams through four loops that corre- 
spond to dual conformal integrals in four dimensions with explicit numerator 
factors of // 2 . (There are no examples of this type below three loops). In the 
bottom row we have isolated three degenerate diagrams which are constrained 
by dual conformal invariance to be equal to pure numbers (independent of s, 
t and // 2 ). 

Evaluation of Dual Conformal Integrals 

In this section we describe the evaluation of dual conformal integrals. Even though 
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d4 



e3 



h 



Fig. 9: Type IV: All remaining dual conformal diagrams. All of the corre- 
sponding off-shell integrals diverge in four dimensions. 

dual conformal integrals are finite in four dimensions we evaluate them by first dimension- 
ally regulating the integral to D = 4 — 2e dimensions, and then analytically continuing e 
to zero using for example algorithms described in [p70| , p7T]JTgf| . For a true dual conformal 



invariant integral the result of this analytic continuation will be an integral that is finite 
in four dimensions, so that we can then freely set e = 0. However the type II and type IV 
diagrams shown in Figs. 6 and 8 turn out to not be finite in four dimensions (as can be 
verified either by direct calculation, or by applying the argument used in [^9| to identify 
divergences). This leaves (1,1,4,17) integrals to be evaluated at (1,2,3,4) loops respectively. 

4-1. Previously known integrals 

Here we briefly review the (1,1,2,5) off-shell dual conformal integrals that have already 
been evaluated in the literature. 

The first class of integrals that have been evaluated off-shell are those corresponding 
to the ladder diagrams JW, J( 2 ), X^ a = 1^ and jW a = JW. In fact an explicit formula 
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for the off-shell L-loop ladder diagram was given in the remarkable paper [HI . The function 
$>( L '(x, x) in that paper corresponds precisely to our conventions in defining the ladder 
diagrams X^ (including the appropriate conformal numerator factors), so we copy here 
their result 



x 



y/l-4x 



Li 2 L (-!/)+ Yl WW ,wr n ; c(fc+oio g 



2L-l-k 



k,l=0 
k-\-l even 



fc!Z!(L-Jfe)!(L-0! 



where 



y 

(4.1) 

, (4.2) 

1 - 2x + VI - 4x 

The second class of integrals that have been evaluated off-shell are those that can be 
proven equal to X^ using the 'magic identities' of [E3|. There it was shown that 



2x 



and 



j(3)o _ j(3)6 = j(3) 



j(4)a = j(4)b = j(4)c = j(4)d = j(4)e = j(4) _ 



(4.3) 



(4.4) 



These identities appear to be highly nontrivial and are only valid for the off-shell integrals 
in four dimensions; certainly no hint of any relation between these integrals is apparent 
when they are taken on-shell and evaluated in 4 — 2e dimensions as in 

Jill- Moreover 



m 



23| the relations (|4.3|) and (|4.4|) were given a simple diagrammatic interpretation which 



can be utilized to systematically identify equalities between certain integrals at any number 
of loops. 



4-2. New integrals 

As indicated above we evaluate these integrals by starting with Mellin- Barnes repre- 
sentations in 4 — 2e dimensions and then analytically continuing e to zero. A very useful 



Mathematica code which automates this process has been written by Czakon |T9|]. We 
found however that this implementation was too slow to handle some of the new integrals 
in a reasonable amount of time so we implemented the algorithm in a C program instead. 
The most difficult off-shell integral we have evaluated, X^ 4 ^ 2 , starts off in 4 — 2e dimensions 
as a 24-fold Mellin- Barnes representation (far more complicated than any of the on-shell 
four-loop integrals considered in [[24 , 25] , which require at most 14- fold representations), yet 
the analytic continuation to 4 dimensions takes only a fraction of a second in C. In what 



12 



follows we display Mellin-Barnes representations for the various integrals in 4 dimensions, 
after the analytic continuation has been performed and e has been set to zero. 

The most surprising aspect of the formulas given below is that we are able to write 
each integral in terms of just a single Mellin-Barnes integral in four dimensions. This 
stands in stark contrast to dimensionally regulated on-shell integrals, for which the analytic 
continuation towards e = can generate (at four loops, for example) thousands or even tens 
of thousands of terms. For the off-shell integrals studied here something rather amazing 
happens: the analytic continuation still produces thousands of terms (or more), but for 
each off-shell integral it turns out that only one of the resulting terms is non-vanishing at 
e = 0, leaving in each case only a single Mellin-Barnes integral in four dimensions. 

This surprising result is not automatic but depends on a number of factors, including 
the choice of initial Mellin-Barnes representation, the choice of integration contour for the 
Mellin-Barnes variables Zi, and some details of the how the analytic continuation is carried 
out. All of these steps involve highly non- unique choices, and by making different choices it 
is easy to end up with more than one term that is finite in four dimensions. However in such 
cases it is always possible to 'reassemble' the finite terms into the one-term representations 
shown here by shifting the contours of the remaining integration variables. 

For the new 3-loop integrals, both shown in Fig. 8, we find the Mellin-Barnes repre- 
sentationsi 



J(3)c = - / ttt-t? xZ2 r(-zi)r(-z 2 )r(z 2 + i)r(-z 2 + z 3 + i)r( Zl + z 3 -z 4 + i) 

T(-z 4 ) 2 T(z 4 - z 3 )T{-z 2 - z b ) 2 T{z 1 - z 4 - z 5 + 1) 
T(zi - z 2 + z 3 - z 4 - z 5 + l)r(-2i + z 4 + z 5 )T(z 2 + z 4 + z 5 + 1) 
T{-z x + z 2 + z 4 + z 5 )T(—zi - z 3 + z 4 + z 5 - l)/(r(l - Zi) 
T(-z 2 - z 5 + l)r(zi - z 2 + z 3 - z 4 - z 5 + 2)r(— z\ - z 2 + z 4 + z 5 ) 

r(— z\ + z 2 + z 4 + z 5 + 1)) 



All formulas in this section are valid when the integration contours for the Zi are chosen to 
be straight lines parallel to the imaginary axis and such that the arguments of all T functions in 
the numerator of the integrand have positive real part. 
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and 

J(3)d = / iP^ r(-zi)r(z! + i)t(-z 2 )t( Zi + z 2 -zs + i)r(-z 3 ) 2 r(^ - *0 

r(^ 2 - 2 3 - z 4 + i)r(zi + 2 2 - z 3 - z 4 + i)r(-z 4 ) 2 r(z 3 + z A + i) ^ 6 ^ 

r(-z 2 + z 3 + Z4)T(-zi -Z2+Z3 + z 4 )/(F(l - zx)T(l - z 3 )T(l - z 4 ) 
T(zi + z 2 - z 3 - z 4 + 2)r(-z 2 + z 3 + z 4 + 1)). 

As expected, T^ d is x-independent because the corresponding diagram is degenerate. 
Upon evaluating ( |4.6|) numerically (using CUBA we find 

j(3)d w 20.73855510 (4.7) 



with a reported estimated numerical uncertainty smaller than the last digit shown. 

Finally we have 12 off-shell four- loop integrals left to evaluate, corresponding to the 9 
diagrams shown in Fig. 8, along with three of the diagrams 

(J (4)/ 5 J(4)d2 and j(4)/2 )from 

Fig. 6. It turns out that 4 of these 12 integrals (J (4)/ , J (4)/ ', j( 4 ) e2 ' and j( 4 ) c ') are signif- 
icantly more difficult than the rest because they apparently require analytic continuation 
not only in e but also in a second parameter v parameterizing the power of the numerator 
factors. (That is, the integrals initially converge only for v < 1 and must be analytically 
continued to v = 1.) We postpone the study of these more complicated integrals to future 
work. 

In analyzing the remaining 8 off-shell four-loop integrals we have found two new 'magic 
identities', 

j(4)e2 _ j(4)W ? j(4)cl _ _j(4)d2_ /^gX 

We established these results directly by deriving Mellin-Barnes representations for these 
integrals and showing that they can be related to each other under a suitable change of 
integration variables. It would certainly be interesting to understand the origin of the 
relations (|4.S| ) and to see whether the insight gained thereby can be used to relate various 
dual conformal integrals at higher loops to each other. 
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Mellin-Barnes representations for the 8 off-shell four-loop integrals are: 

j(4)ci = _ j( 4 )d 2 = _ r xZ2 T{ _ Zi _ 1)r(zi + 2 )r(- 22 )r(^ 2 + 1) 

J (2m)' 

r(-zi - z 3 - i)r(-z 3 )r(z 3 - * 2 )r(zi - * 2 + z 3 + i)r(-z 4 ) 

r(zi - ^ 2 + ^3 + ^ + 2)r(-z 6 ) 2 r(^ 4 + 25-^7 + 1) 

r(-zi + z 2 + z 4 - zq - z 7 )r(z 4 + z 5 - z 6 - z 7 + i)r(-z 7 ) 2 r(z 7 - z 5 ) ( 4-9 ) 

Y(z & + z 7 + l)T(-z 4 + z 6 + z 7 )T(-z 3 - z 4 - z 5 + z 6 + z 7 - 1)/ 

(r(- Zl -Z2- i)r(i - z 3 )r(zi -z 2 + z 3 + 2)r(i - z 6 )r(i - z 7 ) 
r(z 4 + 25 - 2 6 - ^ 7 + 2)r(-z 4 + z 6 + 2 7 + 1)) 



x (4)/2 = | ^^^r(-z 1 -z 1 o) 2 r(-z 2 )r(-^)r(^ 1 +z 1 o + ^ + i) 2 r(-^) 
r(-z 5 )r(-z 6 )r(z 10 + z 2 + z q )y(-z 7 ) 2 y(z 4 + z 7 + 1) 2 

r(— Z\ — Z4 — Z5 — Z 8 — 2)(—Zi — Zio — Z 2 — Z 4 — Zq — Z 7 — Z 8 — 2) 
r(-Z 8 )r(zi + ^ 8 + 2)(Z 2 + -24 + ^5 + ^7 + Z 8 + 2) 

r(-Z 2 - Z 3 - Z 4 - Z 5 - Zg - 2)r(-Zi ~ Z 2 - Z 3 - Zq - Zg)T(-Zg) 
F(z 1 + Z W + Z 2 + Z 3 + Zq + Z 9 + 2) 

T(z 10 + Z 2 + Z 3 + Z 4 + Z 5 + Zq + Zg + 2)T(-Z 1 + Z 3 - Z 8 + Zg) 

T(z 2 + Z 4 + Z 5 + Zq + z 8 + z 9 + 2)/(T(z 1 + z 10 + z 3 + 2)T(-z 4 - z 5 ) 
F(-zi - zio - z 7 )Y(z 4 + z 7 + 2)T(-zi - zio - z 2 - zq - z 8 ) 

T(-Z 3 - Zg)T(-Zi + Zio + Z 2 + Z 3 + Zq - Z 8 + Zg) 
F(zi + Zio + Z 2 + Z 3 + Z 4 + Z 5 + Zq + Z 8 + Zg + 4)) 



(4.10) 



J(4)e6 = " / J^r x ^ r( ~- z ^- z ^ r ^ + i) 2 r(-^)r(^ 3 + 1) 

T(zi + z 3 -z 4 + l)T(-z 4 ) 2 T(z 4 - z 3 )T( Zl -z 4 -z 5 + l) 

r(zi + z 3 - z 4 - z 5 + i)r(-z 5 ) 2 r(z 4 + z 5 + i)v(-zi + z 4 + z 5 ) 

T(-zi -z 3 + z 4 + z 5 )/(T(-2z 2 )T(l - z 3 )T(l - z 4 )T(l - z 5 ) 
T(zi + z 3 -z 4 -z 5 + 2)Y{-zi + z 4 + z 5 + 1)) 
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X (4)e2 = X ( 4 )6i = _ J _dj_ xZ3 r(_ Zl )r(-z3)r(^ + i)r(z 3 - Z 2 )T(-Z 4 ) 2 
r(z 2 -z 3 + z 4 + i) 2 r(-z 5 ) 2 r(z 1 + z 2 -z 5 -z 6 + 2)r(-z 6 ) 2 

F(z 5 + Z 6 + l)r(-Zl - Z 2 + ^5 + z G - i)r(-2i - Z 2 + Z 3 + Z 5 + Z 6 
r(-z 4 + z 5 - z r )T(-zi - z 2 - z± + z 5 + z 6 - z 7 - l)r(-z 7 ) 
r(— Z 3 + Z 4 + z 7 )T(zi + z 2 - z 3 + z A - z 5 + z 7 + 2) 
r(zi - Z 5 - Z 6 + Zj + 1)/ (T(z 2 - Z 3 + Z± - Z 5 + 2)r(-24 - Z 6 + l) 
r(-Zl - Z 2 - ^3 + ^5 + Z 6 - 1)T(-Zi - Z 2 +Z 3 + Z 5 + Zq) 

T(-z 4 - z 7 + l)T(zi + z 2 - z 3 + Z 4 - z 5 - z 6 + z 7 + 2)) 



jim = J( 4 )e i = f r(-z 1 )r(-z 2 )r(-z 3 )r(z 3 + i)r(-^ - z 2 - z 4 - 2) 

7 (27TZj D 

r(-zi - 2 2 - ^3 - Z4 - 2)r(-z 4 ) 

T(zi + z 2 + z 4 + 3)T(zi + z 2 + z 3 + z 4 + S)T(zi + z 3 - z 5 + l)T(—z 5 ) 

T(z 5 - z 3 )T(z 2 + z A + z 5 + 2)T(-z 4 - z 5 - z 6 - l)T(-z 6 ) 2 ( 4 - 13 ) 

r(z 4 + z 6 + i) 2 r(z 2 + z A + z 5 + z 6 + 2)/(r(i - z 3 ) 

T{-z\ — z% — z±- l)T{z\ + z 2 + z 3 + z 4 + 4)r(l - 25 ) 

r(z 2 + z 4 + z 5 + 3)r(i - 2 6 )r(2 4 + z 6 + 2)) 

We do not consider the equality of the two degenerate integrals X^ b2 = X^ el to be 
a 'magic' identity since it is easily seen to be a trivial consequence of dual conformal 
invariance. Evaluating them numerically we find 

j(4)62 =J (4)el = 7 o.5 9) (4.14) 

again with a reported estimated numerical uncertainty smaller than the last digit shown. 

It would certainly be interesting to obtain fully explicit analytic results for these new 
integrals. Although this might seem to be a formidable challenge, the fact that it has been 



possible for the ladder diagrams (|4.1| ) suggests that there is hope. 
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4-3. Infrared singularity structure 

Finally, it is clearly of interest to isolate the infrared singularities of the various inte- 
grals. For the previously known integrals reviewed in subsection 3.1 we expand (|4.1|) for 
small x, finding 

Z« = log 2 £ + 0(1) 
1 7T 2 

J (2) = - log 4 x + — log 2 X + 0(1), 

1 r 2 7 4 (4.15) 

J(3) = _L l og 6 x + lQg 4 X + J^_ ! 2 

36 6 36 6 36 5 w ' 

^41 1 , 8 77T 2 , 6 497T 4 , 4 31tT 6 , o 

2W = _ log 8 * + — log 6 X + — log" x + — log 2 x- + 0(1). 

For the new integrals evaluated in this paper we obtain the small x expansion directly from 
the Mellin-Barnes representations given in section 3.2 by writing each one in the form 

^Fiy), (4.16) 

shifting the y contour of integration to the left until it sits directly on the imaginary axis 
(picking up terms along the way from any poles crossed), expanding the resulting integrand 
around y = and then using the fact that the coefficient of the l/y k singularity at y = 
corresponds in x space to the coefficient of the 1 fc / log x singularity at x = 0. In this 
manner we find 



C(3) 



4 



j(3)c = M_l log 3 x _ log 2 x + 14 _ 3 2388625 log x + 0(1) (4.17) 
at three loops and 

J( 4) d 2 = _ j( 4)cl = _Cg) log 5 x + 7^ lQg 4 x _ 6 . 75193310 lQg 3 x 

+ 15.45727322 log 2 x - 41.26913 log x + 0(1), 

^<4)f2 1 , 8 77T 2 6 1497T 4 , 4 

X(^ = _, og 8 x + — log 6 x + — log 1 * 

+ 64.34694867 log 2 x + 0(1), ( 4 - 18 ) 
j(4) e6 = _20.73855510 log 2 x + 0(1), 

J(4)e2 = X (4)6i = _ JL_ i og 4 x + i. 7 2821293 log 3 x 
720 5 5 

- 12.84395616 log 2 x + 52.34900 log x + 0(1) 
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at four loops, where some coefficients have only been evaluated numerically with an esti- 
mated uncertainly smaller than the last digit shown. Interestingly, the coefficient of log 2 x 
m 

J( 4)e6 

appears to be precisely (minus) the value of X^ d shown in (|4.71 ). Perhaps this 
can be traced to the diagrammatic relation that is evident in Fig. 8: X^ d appears in the 
'upper diagonal' of j( 4 ) e6 . 

5. Summary 

We have classified all four-point dual conformal Feynman diagrams through four 
loops. In addition to the previously known (1,1,2,8) integrals (Fig. 6) that contribute 
to the dimensionally-regulated on-shell amplitude respectively at (1,2,3,4) loops, we find 
(0, 0, 2, 9) new dual conformal integrals (Fig. 8) that vanish on-shell in D = 4 — 2e but not 
off-shell in D = 4. There are also (0, 0, 0, 11) dual conformal diagrams (Figs. 7 and 9) that 
diverge in four dimensions even when taken off-shell and therefore do not give rise to true 
dual conformal integrals. 

Next we addressed the problem of evaluating new off-shell integrals in four dimen- 
sions. Of the total number (1, 1, 4, 17) of such integrals, explicit results for (1, 1, 2, 5) have 
appeared previously in @,^3[. We find Mellin-Barnes representations for an additional 
(0,0,2,8) integrals, including two pairs related by new 'magic identities', and evaluate 
their infrared singularity structure explicitly. Evaluation of the remaining (0, 0, 0, 4) inte- 
grals j( 4 ^, j( 4 ^', j( 4 ) e2 ', and j( 4 ) c ' is left for future work. 
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